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Abstract 

We study the normalization of perturbative QCD corrections to the inclu- 
sive B Xg^ decay. We propose to set the renormalization scale using 
the Brodsky-Lepage-Mackenzie (BLM) method. In the proposed method the 
scale is determined by absorbing the vacuum polarization correction from 
light fermions to renormalization scale but not including the anomalous di- 
mensions. The BLM scale depends in general on the renormalization scheme 
and the factorization scale. We find that the BLM scale is insensitive to the 
factorization scale. In the heavy-quark potential scheme, we find that the 
BLM scale is Hblm ~ (0.315 — 0.334)mfc when the factorization scale varies 

from mi)/2 to 2mb. 
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The rare decay process b ^ s'y has been recently observed at CLEO II |T[] at a branching 
ratio 2.32(±0.57 ± 0.35) x 10~^ (statistical and systematic errors). The process serves as 
an excellent probe to the new physics beyond the Standard Model such as charged Higgs 
0, supersymmetry 0] or anomalous WW'j coupling [Q. At leading order the process is 
described by an electromagnetic penguin diagram. It is known that there are large QCD 
corrections to the penguin diagram. A complete renormalization group- improved perturba- 
tive calculation is now available to the leading-logarithm approximation. It is found that 
the QCD corrections increase the rate by a factor two to three. In the Standard Model the 
branching ratio to leading-logarithmic accuracy is 2.8(±0.8) x 10^^, where the error is dom- 
inated by the uncertainty in the QCD renormalization scale varied over mi,/2 < fx < 2mb 
I^J^. Inclusion of the next-to- leading logarithmic QCD corrections will reduce the theoreti- 
cal error. However, a complete calculation is still not achieved. With partial next-to-leading 
corrections known so far, the branching ratio falls to ~ 1.9 x 10~^. 

In view of the current status of the theoretical uncertainties as summarized above it is 
of interest to normalize the leading order QCD corrections and estimate theoretical errors 
thereof. In this paper we address this issue and determine the renormalization scale accord- 
ing to the Brodsky-Lepage-Mackenzie (BLM) method |^. So far applications of the ELM 
method have been mostly restricted to processes involving (partially) conserved currents. In 
this case QCD correction may be calculated by a perturbative expansion at a fixed order. 
For example, BLM scale setting has been established for the inclusive semileptonic B decay 
B XqCV H. In Ref. || the BLM scale in the WS scheme is found ^.blm ~ 0.07 
(0.12 when running MS b-quark mass is used), indicating a significant QCD cor- 
rection from higher order terms. The QCD corrections for the radiative B decay is more 
complicated since the process involves large logarithms arising from QCD effects between 
the top-quark or W-gauge boson mass and the bottom quark mass scales. The effective 
Hamiltonian approach takes care of these large logarithms by introducing a factorization 
scale and separates the long- and the short-distance QCD effects. 

In situations involving large logarithms a physical observable P such as the total decay 
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rate or the branching ratio may be written schematically as 



TT 



+ ••• 



(1) 



Here is a factorization scale dividing long- and short-distance physics and hr is a renor- 
malization scale for the QCD corrections. All the large logarithms of {m.w/lJ') are resummed 
into M{iJ,). QCD corrections to the matrix element at n, Hr ^ ^qcd are given by a per- 
turbative expansion. The coefficients Tq, ri and r2 depend on both the factorization and the 
renormahzation scales. Of course the physical observable P calculated to all orders should 
be independent of the two scales. However, at any finite order, perturbative evaluation of P 
do depend on the factorization and the renormahzation scales. In practice the factorization 
scale is set equal to the renormahzation scale, and then the theoretical error is estimated. 
To be precise these two scales are distinct, and there is no reason why the two scales should 
be set equal. For example, in the conformally invariant QCD, a finite-order perturbatiion 
expansion does not depend on any renormahzation scale but still depends explicitly on the 
factorization scale. 

We are interested in calculating perturbative QCD corrections and setting the renormal- 
ization scale /ir for a given renormahzation scheme and the factorization scale /i. Typically 
there are logarithms of {/ir/ ii) present from two sources. One is the logarithms that can 
be resummed to the running coupling constants. Another is the logarithms that may be 
resummed into the anomalous dimensions of operators involved. In order to set the renor- 
mahzation scale we propose to follow the original spirit of BLM closely. The idea is to set /jLr 
by absorbing UfO^ contributions of the vacuum polarization to the renormahzation scale, 
but not including the contributions of the anomalous dimensions. Then the generalized 
BLM scale is set by the fixed-order perturbative correction part inside the bracket in Eq.(l): 
fJ^BLM = fJ'R exp(3r2/ri). If the content of the theory does not change between /x and iir, 
the normalized physical quantity P may be given by 
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P = M(/iBLM) IroifiBLM) + riiflBLAl) '^''^^^^^^ ] (2) 



TT 



in which the anomalous dimension effects are talcen into account by running from M{fi) to 
M{^blm) using the renormahzation group equation. Aside from setting the BLM scale it 
may be also argued that the fixed order coefficients ri,r2 are the most significant parts in 
the next-to-leading order contribution 

We emphasize, however, that our proposal is by no means unique or superior to other 
plausible generalizations of the BLM method. For example, we may set the BLM scale 
including contributions from the next-to-leading order anomalous dimensions as well as ef- 
fective coupling constant through the vacuum polarization. This is complicated and in order 
to be consistent, we need anomalous dimensions at next-to- leading order. Furthermore at 
present a complete next-to-leading order calculation of the coefficient functions and anoma- 
lous dimensions is not available for the radiative B decay. On the other hand our method 
provides a physically motivated generalization of the original BLM method and the final 
result shows a reasonable behavior. 

For the radiative B decay, the effective Hamiltonian below Mw scale is given by ||10[ 



^eff = -Vt,VZ % jz aif^) O.ifi), (3) 

where /i is the factorization scale separating the coefficient functions and the operators. For 
the partonic decays 6 — S7 and b — > s'yg, the dominant contributions in the radiative decay 



arise from the operators |T0 



O2 = {cLal^h Loci's LplfjiCLp) 



e 

16^ 
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^7 = T7r^F^,-s^(j^'''{m,PR + m,Pi)&o 



08 = -^G^^-s^a^^T^^{m,Pa + m,Pi)6^ (4) 

where are color indices, Pr_l = (1 ± 75)/2 and Ffj,i^,G^^, are photon and gluon field 
strengths respectively. Contributions of the other operators enter through operator mixing 



at next-to-leading order and are numerically negligible |TT|. Because of color structure, Oi 
does not contribute to the real gluon emission to the order we calculate. In addition Cs{mb) 
is small compared to the others. Hence we neglect Oi, Og contributions in what follows. 
The Wilson coefficients at leading order are given by |jlO 



^-6/23 _^ ^12/23 



C2imw) 

58 



29 
189 



^28/23 



(5) 



with Tj = as{fi) / asi'iTT'w) and 



C2{mw) = 1, 



X 



Cjimw) = ^^fzY^^^^^^^ - 2) In a; 
- (x- l)(8x=^ + 5x-7)], 

with X = mf/mly. 

Keeping only O2 and O7 the total b ^ decay rate can be written as 



(6) 



r?? + r22 + r 
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where r22(r77) denotes the decay rate obtained from the matrix elements squared for O2 
(Oj) and is the decay rate from the interference terms between O2 and O^. Including 
both virtual correction and real gluon emission r77 is given by [|l^ 



77 



^ 2a,(/ifl) 27r^ 8 

Stt ^ 3 3^ " 



(8) 



where Fq = |At pQ;G|.m^/(327r^) and Xt = VtbV^l. Here we neglect the s quark mass and the 
ellipsis denotes higher order contributions in a^- 

To set the generalized BLM scale according to our proposal, it is necessary to calculate 
the QCD corrections of order a^, which are proportional to the number n/ of the light 
quark flavors. Smith and Voloshin [|12| have recently shown that the nj-dependent part of 
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the order contribution may be written in terms of the one-loop corrections evaluated with 
a fictitious gluon mass A 



X 



r«(A) 



mt 



r(i)(o) 



A2 



(9) 



where b = 11 — 2nf/3 and r(^^(A) is the order as contribution to the decay rate computed 
with a finite gluon mass A and a^g^\mb) is the QCD coupling obtained in the heavy-quark 
potential scheme. 

In order to obtain Tji we calculate the virtual and bremsstrahlung corrections with a 
finite gluon mass coming from O7. The virtual massive gluon correction is obtained by 
calculating diagrams corresponding to the dressing of the vertex, the initial-state b quark 
and the final-state s quark line with a massive gluon. The bremsstrahlung correction is 
calculated from the emission of a massive gluon from either of the external b, s quarks. Sum 
of the two contributions are infrared finite. The sum is calculated analytically as 



a. 



X 



3n 
16x 



rl6 4 2 
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-X H — X —ox mx 
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x(3x^ - lOx - 20) _i (1 - a;) Jx(4 - x) 
H , tan 
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where x is the dimensionless gluon to b quark mass ratio, x = X^/ml and Fq is the tree level 
decay width. 

After carrying out the integration of Eq. (|10|) over x according to Eq. Tjj at order 
al is given by 

4a^{mb),n'' 4. 



T77 = {Cr{fi)yTo 



1 - 



3n 



x(l + ^^np)[2.733] 



(11) 



We can apply the same method for the calculation of and r27. There are no virtual 
corrections at order a^. The double differential decay rates at order with a finite gluon 
mass are given by 
cir22 2as 



dXqdx^ 27 TC 



ro(C2(At))^ X [2{1-Af - X - x^) + xxg) 



— 2(1 — y4^)|x^(l — X — X^) + Xg{l + X — Xg) + Xq{l — X — X 
+ (1 + y4)^|x^(l + X — Xg) + ^ —{X^{1 — X — X^) + Xg{l + X — Xg 

f . .2 8(1 4(1 - A) V ^ 
+ |-1 -2A + + '— - 

A. lAj Q Y -'- i-t. Ju q J 

X |-X-y(l + X - Xg) + ^''' (X^(l - X - X^) + X<;(1 + X - Xg))} , (l2) 

where Xj = 2Ei/mb {i = q^ljQ are s quark, photon, gluon in the final state) satisfying 
Xg + x^ + Xg = 2 and x = A^/m^. And 

A= In^^. (13) 

Similarly 

_|ji_._^r„c,(,)c,(,) 



— — f2(l + X - x„)(l - X - x„) + x^(l - X - x^) + + ^ — ^ 

— Xq{l — X — Xq) — Xg(l + X — Xg) — 4x(l + X — Xg) 



+ A[{^Xg — 4X^)(1 + X — Xg) + (4(1 + X — Xg) + Xg)(l — X — Xg) " X^(l " X " X^ 

+ ^(1 + a; _ a; )(i _ a; _ :r ) _ 2xx^£+x^Xg) x . 

1 1 

I ^ ^ gr ^ 1 ij/ ^ g ^ I OC-^ ^ 1 ^ X ■-y ^ ^ p ^ 1 I X X g ^ | /y 

1 X'y 1 



+ a{—2xx^—^ H ^(3^0(1 + X — Xn) + x^(l — X — x-v) — Xo(l — X — x„)) 

— 4x^(1 + X — Xg) — X5(l — X — Xg) — x^(l — X — x-^) + Xg(l + r — Xg)^ I . (14) 

These are infrared finite as x — >■ and the phase-space integration and the integration 
over the gluon mass yield finite results. We can get the single differential decay rates by 
integrating Eqs. ([T2|) and ( p!^ over, say Xg, analytically and the remainingg calculation is 
done numerically. Using the Eq. (§) and after the integration, we find 



^22 
^27 



2-L 0- 
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[2.486] + ■ ■ 
[6.648] + 



(15) 
(16) 



Summing Eqs. (PD , ([TsD, and ( |I6D the total decay rate is given by 

r(/i) = T77{n) + r22(/i) + r27(/i) 

r / 

= ro(C7(/i))^ 1 + — ( -7.83 - 0.222k + O.USk' 



/ ha, 21.5 + 1.48k - 0.368^2 
(1 + ^ ^., . .... ^^^) + 



(171 

47r 7.84 + 0.222k -0.148^2^ ' J' ^ ^ 

where k = C2/C7 is the ratio of the Wilson coefficient functions. In Eq. ([TtD the terms 
proportional to can be absorbed into the order as term according to the BLM prescription. 
Recall that there are two independent scales in the decay rate. One is the factorization scale 
/i at which the Wilson coefficient functions and the local operators are factorized and the 
other is the renormalization scale hr at which the theory is renormalized, and the two scales 
are independent. 

As we have mentioned, we propose to set the BLM scale only from the fixed-order 
perturbation series in asifiR). This method extends the original idea of Brodsky, Lepage 
and Mackenzie [^ closely in the sense that the BLM scale is determined by absorbing the 
vacuum polarizations of light fermions into the coupling constant asifJ^BLAi)- In this case 
the coefficients in front of the powers of ^^(/ii?) are functions of the factorization scale 
through the anomalous dimension of the operators. Therefore the BLM scale depends on 
the factorization scale. Noting that in the heavy quark potential scheme 



Air ml 

the BLM scale from Eq. (|1^ is given by 

1 21.4 + 1.48k - 0.368k2 

AisLM = mb exp 



(18) 



(19) 



2 7.83 + 0.222k -0.148k2- 

Here k is a function of the factorization scale /i. In Fig. 1, the BLM scale hblm is plotted 
as a function of the factorization scale /x. For numerical evaluation we have used = 180 
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GeV, mw = 80.1 GeV, mb = 5 GeV and A = 225 MeV for nj = 5. Our result in Fig. 1 
shows that the BLM scale is extremely insensitive to the factorization scale. As we vary 
the factorization scale fi from mb/2 to 2mb, the BLM scale changes only over 0.315 to 
0.334 mfe. This is gratifying since the insensitivity of the factorization scale implies that we 
can reduce part of the theoretical error of the decay rate from the choice of the factorization 
scale. In our case the renormalization scale is almost constant and it yields the normalized 
decay rate which does not vary widely when the factorization scale changes from mb/2 to 
rub- 

The total decay rate is then normalized by scaling down the leading order coefficient 
functions from /i to /Ublm using the leading logarithmic coefficient function. According to 
the original idea of BLM, the normalized rate may be interpreted as the mean total decay 
rate averaged over the gluon virtuality. In order to compare this normalized total decay rate 
with other decay rates, we have plotted three decay rates as a function of the factorization 
scale in Fig. 2. The dotted line represents the a°-order result from Oj only. The dashed line 
is the as order result in Eq. (^) with the factorization scale /i and the solid line is the 
result with /i = ^BLMif^)- As shown in Fig. 2, the result with the BLM scale is almost flat in 
the range while other results vary. Note that in plotting Fig. 2, we have used the coefficient 
functions only to the leading logarithm order. If the complete next-to-leading logarithmic 
results for the coefficient functions and the anomalous dimensions are available in the future, 
it will be interesting to compare the next-to-leading log result of the decay rate with our 
BLM result which will also be improved by using the coefficient functions at next-to-leading 
orer. This will eventually test whether our proposal of setting the renormalization scale is 
sensible. 

We can also ask questions like how the BLM scales and the decay rate can be expressed 
in other schemes. Since there exists a definite relationship of ^blm and physical observables 
among various schemes, our result can be converted without any ambiguity to other schemes 



such as the MS scheme and can also be extended to include the running mass effects of the 
h quark, [l^ 
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FIGURES 

FIG. 1. The BLM scale HBLM/nT-b for the total decay rate F as a function of n/mi,. 



FIG. 2. The total decay rates F/Fq as a function of n/mi,. 
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